Effective Quantum Dynamics of Interacting Systems with Inhomogeneous Coupling 
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We study the quantum dynamics of a single mode/particle interacting inhomogeneously with a 
large number of particles and introduce an effective approach to find the accessible Hilbert space 
where the dynamics takes place. Two relevant examples are given: the inhomogeneous Tavis- 
Cummings model (e.g., N atomic qubits coupled to a single cavity mode, or to a motional mode in 
trapped ions) and the inhomogeneous coupling of an electron spin to N nuclear spins in a quantum 
dot. 

PACS numbers: 42.50.Fx, 42.50.Vk, 73. 21. La 



I. INTRODUCTION 

In quantum optics and condensed matter physics, a 
great effort has been oriented to the study of single quan- 
tum systems, or a few of them, and their experimental 
coherent control. This is the case of cavity QED 
trapped ions 0, quantum dots 0], circuit QED 
among others. On the other hand, many-particle physics 
describes a variety of systems where collective effects may 
play a central role. 

The interaction of a single atom with a quantized elec- 
tromagnetic mode, described by the Jaynes-Cummings 
(JC) model has played a central role in quantum op- 
tics and other related physical systems. Among its many 
fundamental predictions, we can mention vacuum-field 
Rabi Oscillations, collapses and revivals of atomic popu- 
lations, and squeezing of the radiation field [6j]. Atomic 
cloud physics can be conveniently described by the Dickc 
model [7], when considering the interaction of atoms 
with light in free space, or by the Tavis-Cummings (TC) 
model [8] , when the coupling takes place inside a cavity. 
In most treatments and applications of the TC model, a 
constant coupling between the atoms and the radiation 
field is assumed, simplification that greatly reduces the 
complexity of an analytical description Q EH EH, E3j- 
Even in the case of many atoms, the homogeneous- 
coupling case profits from the ST/ (2) group structure of 
the collective operators, allowing at least numerical solu- 
tions. However, the situation is drastically different when 
we consider the more realistic case of inhomogeneous cou- 
pling. In this case, there is no such straightforward way 
to access the Hilbert space, because all angular momenta 
representations arc mixed through the dynamical evolu- 
tion of the system and no analytical approach is known. 

On the other hand, electron and nuclear spin dynamics 
in semiconductor nanostructures has become of central 
interest in last years 0, EH, an d several techniques for 
quantum information processing have been proposed [hH . 
Ha . EE E3, El, El]- For example, ensembles of polarized 
nuclear spins in quantum dots 12011 and quantum Hall 
semiconductor nanostructures [2l| have been proposed 



as a long-lived quantum memory for the electron spin. 
Recently, polarization procedures for the nuclear spins 
in a q uantum dot have been studied theoretically (22l . 
l23l. |24|| . and first experiments have achieved considerable 
dynamical nuclear polarization [25| . However, difficulties 
in producing a fully polarized nuclear state suggest the 
necessity of analyzing the effect of imperfect polarization. 
This treatment faces the same theoretical challenges as 
the study of the inhomogeneous TC model, because of 
the unavoidable inhomogeneous coupling of electron and 
nuclear spins in quantum dots. 

In this work, we develop an effective approach to study 
the dynamics of a single mode /particle coupled inhomo- 
geneously to a large number of particles. In particular, 
we analyze in Sec. |H] the general interaction of a col- 
lection of two-level atoms with a quantized field mode, 
that is, the inhomogeneous TC model. We identify a 
subset of the Hilbert space that is relevant for the dy- 
namics and show that we are able to reproduce accurate 
results with a small number of states. In Sec. IHII we con- 
sider the hyperfine interaction of a single excess electron 
spin confined in a quantum dot with quasi-polarized nu- 
clear spins. Using the formalism developed in section HT1 
we study the transfer of quantum information between 
those systems in the presence of a single excitation in 
the nuclear spin system prior to the write-in procedure. 
In Sec. IV, we present our concluding remarks. 



II. INHOMOGENEOUS TAVIS-CUMMINGS 
MODEL 

The Hamiltonian describing the inhomogeneous cou- 
pling of N atoms with a quantized single field mode, in 
the interaction picture and under resonant conditions, 
can be written as 
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Here, gi — g(fi) is the (real) inhomogeneous coupling 
strength of atom i at position r^, cK~ {erf) is the lowering 
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(raising) operator of atom i, and a (a^) is the annihila- 
tion (creation) operator of the field mode. We will refer to 
this model as the inhomogeneous Tavis-Cummings (ITC) 
model. In the homogeneous case, gi — g, Vi, one can de- 
fine as angular momentum operator J + oc ffX^iLi^j 
describing transitions between the common eigenvectors 
of J z and J 2 , where J z — X)i=i This is not true for 

the inhomogeneous case, because J + = 5Z i=1 g^t does 
not satisfy anymore the SU (2) algebra. However, we will 
show that this problem can still be addressed with cer- 
tain restrictions. The strategy consists in following the 
Hilbert space that the system will visit along its evolu- 
tion, and implementing a truncation that will depend on 
the window evolution time. 

Let us first consider the initial condition given by 



|*(0))=5> n |n}|0> 



(2) 



where \n) denotes a photonic Fock state and |0) the col- 
lection of N atoms in the ground state. For each product 
state |n)|0), we have a fixed number of excitations, which, 
even in the inhomogeneous case, is conserved by the dy- 
namics determined by Hnc- If the initial state is |0)|O), 
the unitary evolution is trivial 
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Starting from initial state |1)|0), -Hitc produces a non- 
trivial dynamics via the term J + a, so that 
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J+a|l>|0)=|0>$>|li>, 



(4) 



where = o^O) represents a TV-atom state with one 
excitation in the ith atom. We define 
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a normalized state such that 
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The generated state |1)|1) produces now 

N N 
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At the same time 

j-at|i)|i)=v^JV |2)||o). (n) 

On the other hand, the state |0)|2) evolves according to 



j-at|0}|2) = |1) Y,Y,^9i9oK\ 2 i 



= |l)|0i>, 
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with 
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and Cij — 2gigj/N\. The result is a one excitation state 
in the atomic system which is different from |1) defined 
in Eq. ([5]). The resulting state can be expressed as a 
linear combination of |1) and an orthogonal state 
that can be conveniently expressed as 



<^xi^x> - i<ii^i>r 
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yielding 
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Note that the system evolves in a closed two-dimensional 
subspace {|1) |0), |0) |1)} in which all states have the same 
fixed number of excitations. 

If we start the dynamics with the 2-excitation state 
1 2) 1 0) , it evolves through 



J+a|2)|0) = V2\l)^2 9i \li 



i=0 



\/2Ar |l)H) 



(7) 



11)100 = <ii^ >ii)|i)+v<<h<m - i(ii0i>rii)ii P ). 

(15) 

The state |1)|1 P ) is a new orthogonal state of the con- 
sidered 2-excitation Hilbert subspace. Via J~a\ |l)|l p ) 
goes back to |2)|0), the latter being coupled only to the 

state 1 1) 1 1). We observe also that, while J~a) | 1) | l p ) = 0, 
the associated Hilbert subspace grows according to 

J + a\l)\l p ) =a 2 |0)|2)+a 2p |0)|2 p ), (16) 

where we introduced \2 p ) perpendicular to |2). The state 
|2 p ) thus obtained is coupled via J" to states |l p ) and 
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FIG. 1: Schematic representation of the growth of the Hilbert 
space associated with the collective atomic states. 

\lp p ), where |l pp ) is orthogonal to and so on. In 
fact, applying (J^a) n on |n) |0) yields the state 

|0)|n) =(n!) 3 / 2 |0) J2 •■>■ •'/--••//. ....... • (17) 

i\ <...<in 

and the collective atomic states couple following the table 
of Fig. [1] In general, the states |ro)|m) lead to linearly 
independent subspaces of conserved n + m excitations. 

Eventually, all states of the Hilbert space will be visited 
by an arbitrary system evolution but, within a given time 
window, only a restricted portion of the Hilbert space will 
be accessed. However, in the example of Fig. [2] we ob- 
serve that our approximation for the evolution of the col- 
lective atomic ground state is in good agreement with the 
exact result for very long times, reproducing even pop- 
ulation collapses and revivals. Surprisingly, to achieve 
this level of accuracy, it is enough to consider the first 
two rows displayed in Fig. [1] The number of columns 
that has to be considered is determined by the number 
of initial excitations in the system because the Hamilto- 
nian of Eq. |T]) preserves this number. In consequence, 
the low mean photon number, n = 1.8, suggests that the 
initially unpopulated atomic state may reach at most a 
similar mean value. Given that the number of atoms is 
larger, N = 6, it is expected that the normalized cou- 
plings between the few visited atomic states in the first 
row of Fig. [1] and the rows below are quite small. In order 
to confirm this intuition, we show in Fig. [3] how the pop- 
ulation of the second and third rows of Fig. [T] decreases 
as the number of atoms N increases for the same param- 
eters of Fig. [2j Clearly, also increasing inhomogeneity of 
the coupling causes increased leakage into higher order 
rows, see Eq. (fT4")l . The example of Fig. [5] shows a reduc- 
tion of the Hilbert space dimension from 2 6 x 7 = 448, 
where 7 is the dimension of the truncated field space, to 
an effective size of 49. In general, to improve the accu- 
racy, we only need to take into account more rows in the 
calculation. 
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FIG. 2: Evolution of the collective atomic ground state pop- 
ulation, -P|o)> in the case of inhomogeneous coupling gj = 
(/sin [jn/ (N + 1)], j — 1, .., N, with N — 6, an initial atomic 
state |0), and an initial coherent state in the field with 
n = 1.8. Exact evolution (solid line) and approximated evo- 
lution (dotted line). 




FIG. 3: Evolution of the total population of the second and 
third rows, Pr 2 ,r a , for different number N of atoms. 



We have proposed a method to follow the dynamics 
of a system with inhomogeneous coupling in a relatively 
simple manner. It is now interesting to see how strong 
can be the effects of inhomogeneity on the dynamics of 
relevant observables. For example, in Fig.[4j we compare 
the different predictions for the fluctuations of the field 
quadratures. For the homogeneous case, the first time 
scale corresponds to y/Ng, while for the inhomogeneous 
case the analogous would be (J2j g 2 ) 1 ^ 2 ■ In this man- 
ner, to compare both cases, the homogeneous and the 
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FIG. 4: Evolution of the quadrature fluctuations (AXi) — 
(X\) ~ (^i) 2 f° r N = 8, initial atomic ground state, and 
a coherent field with n = 1. Solid line: inhomogeneous 
coupling gj — g sin (N + 1)], where j labels the posi- 
tion of the j-th atom. Dashed line: homogeneous case with 

3 Si, 



g 3 = 1/(V5V)(£ ^ 1/2 



inhomogeneous we consider the homogeneous coupling 
case with coupling g 3 = l/( \/~N)(J2j ffj) 1 ^ 2 - ^ n Fig- HI we 
observe that not only does the inhomogeneous situation 
differ from the homogeneous in the typical timescale of 
the dynamics, but also it shows additional effects, which, 
in the case we consider, are reflected in an increase of the 
quadrature fluctuations. Our method allows the study of 
these effects for particle numbers that are non tractable 
for standard simulation methods. 



III. QUANTUM DOTS 

Another interesting physical system involving inhomo- 
geneous coupling in a natural way is that of a single elec- 
tron confined in a quantum dot surrounded by nuclear 
spins. In presence of a magnetic field Bq along the z-axis 
the associated Hamiltonian, neglecting dipolar interac- 
tions between nuclear spins, can be written as [26l |27| 



jv 



Hqd — 9* HbBqS z + g„/i s 



(i) 



Vhf. 



(18) 



i=l 



Here, the first two terms are the Zeeman energies of 
the electron and nuclear spins, and Vhf = Vd +Vq de- 
scribes the hyperfine contact interaction of a single elec- 
tron spin interacting with N nuclear spins, where Vd = 

Eli ohWSz, and % = EiIiK/2) (i^s+ + ifs-). 

The coefficients on — Avq \tp (fi)\ 2 are the inhomogeneous 
coupling constants, where A is the hyperfine coupling 
constant, vq is the inverse density of nuclei in the mate- 
rial, and ip (r*i) is the envelope function of the electron at 



position Ti. The term Vd induces an effective magnetic 
field for the electron B e g = B + (l/g*^B) YliLi a i^\ 
the so called Overhauser shift. Below, the nuclear Zee- 
man energies will be neglected, as their magnetic moment 
is much smaller than the one for the electron 26]. 

By applying an appropriate field Bq, so that 
g*/iBBes <C Vq, one may allow the spin-exchange term 
to dominate the dynamics [20 ]. In this case, and defining 

A.±,z = J2j a jl± z; t ne Hamiltonian of Eq. (|18[) approxi- 
mately reduces to Vn = (1/2)(A_S'+ + A + S-). However, 
in general, the term Vd = A Z S Z cannot be neglected. 
Even on resonance, where -Bo = — 0-1 9* L l b){A z ), one has 
(A 2 Z ) — (A z ) 2 7^ 0. In conclusion, the Overhauser field felt 
by the electron spin cannot be fully compensated. To in- 
clude this effect we rewrite the zz-part of the interaction 

as A z = aJ z + 5A Z = aj^li + J2i( a i ~ so 
that the on-rcsonancc Hamiltonian becomes 



+ +aS z {j z -{J z ) Q ) + V z , 



(19) 

Here (J 2 )o is the expectation value with respect to the 
fully polarized state, a = EiOti/N is the average cou- 
pling constant (which equals A/N due to the normaliza- 
tion of the electron wave function), and V zz — S Z 5A Z . 
The latter term, i.e. the inhomogeneous contribution to 
the z-term of the interaction, can be treated perturba- 
tively. As we show later, the homogeneous zz-term S Z J Z 
is neglected with respect to the flip-flop term due to a 
small factor 1/^/N. The term V zz is even smaller and 
negligible, because it depends additionally on an inho- 
mogeneity parameter, like the variance of the (smooth) 
electron wave function. 

In this work we consider 7=1/2 nuclear spins in a 
spherical quantum dot, and 



oc e 



(20) 



as the electron envelope function 26] of size r . 



A. Electron and Nuclear Spin Dynamics 

The inhomogeneous coupling of the electron spin to the 
nuclear spins does not allow for a general analytical solu- 
tion of the dynamics. The assumption of a fully polarized 
initial nuclear spin state reduces the dimension of the rel- 
evant Hilbert space from 2 N+1 to JV + 1 and therefore 
the problem is readily analyzed. On the other hand, a 
more general initial state of nuclear spins vastly increases 
the dimension of the Hilbert space that the system will 
visit along the evolution. Motivated by the proposals for 
quantum information storage via the hyperfine interac- 
tion [20| . we will study situations of imperfect, but high, 
polarization. We analyze the case where the dynamics 
involves one excitation in the nuclear spin system, which 
we will refer to as "defect" . We show, with methods sim- 
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ilar to the ones developed in Sec. II, that the relevant 
Hilbert space remains still small. 

The most general pure state describing this situation 
can be written as 



N 



(21) 



with ^ 



N 
3=1 



= 1, and 111 



represents the state 



\H ... ]j ... II) with an inverted nuclear spin at posi- 
tion Tj . In particular we consider the case of a uniformly 
distributed defect Oj = 1/y/N and a localized defect at 
position jo , 



aj =C{T/2f/ (j-j ) 2 + (T/2) 2 



(22) 



where C is a normalization constant and T characterizes 
the width of the distribution. 

Let us consider two possible initial states for the elec- 
tron and the nuclear spins with a single imperfection, 



IT} £ 
IT> £ 



111), 
IT l). 



(23) 
(24) 



In order to solve the dynamical evolution of the compos- 
ite system, we employ the recipe of following the rele- 
vant Hilbert space, as described in Sec. II. We observe 
that, under the action of the Hamiltonian of Eq. (fT9]) , 
the initial state || 1) is coupled to the state |f 0) via 
S + A_ || 1) = 7 IT 0), with |0) = |0)„ = |U ... Il) n and 



7 



(25) 



The state |T 0) will be coupled through the term S-A+ 
to a state with one excitation, but different from |J. 1). 
Similar to the formalism developed in Sec. II, we can 
write S-A+ |T 0) = 7 \l 1) + \l lx), where state. 




a% 1» 



7|1) 



(26) 



is orthogonal to state |1), with (3 = y (No) — 7 2 and 
^0 = J^2 a i- Finally, we observe that S + A- |J, lx) = 
P\l 0). In summary, the exchange terms yield 

S+A-\i 1) = 7 |T0), 
S_i+|T0) = 7 |11)+ J 9|11 X >, 

1 ± ) = /5f |T 0> . (27) 

Furthermore, the zz-interaction yields matrix elements 
(il|a5,(J,-<J,) )|il) = -f, 



(l0\aS z (J z -(J z ) )\l0) = 0, 
(I l±\aS z (j z -(J z ) )\l 1±) : 



Q 

2' 



(28) 



Therefore, for the initial condition of Eq. (|23[) the 
system evolves only in the subspace spanned by 
{|T 1) ,|T 0) ,\l 1 The temporal evolution in this 

subspace can be exactly solved, obtaining (t)) — 
ai(t) |T l) + h(t) |T 0)+Cx(t) \l l x ), where the probabil- 
ity amplitudes are given by 



01 (t) = 1 + 



7 2 exp (iat/A) 
W~8 



8 cos (St/2) 



h(t) 
Cl (t) 



+i— sin(<Si/2) - 5 exp (-iat/A)], 

7 _ 
— exp (iat/A) sin (St/2) , 




7/3 exp (iat/A) 



n 2 s 

-S exp (—iat/A)] 



S cos (St/2) + i- sin (St/2) 



(29) 



Here, fl = y/j 2 + f3 2 = N and S = sJ(a/2) 2 + tt 2 
plays the role of a generalized Rabi frequency for a pro- 
cess with detuning a. The zz-interaction thus leads 
to a small detuning for the exchange of excitation be- 
tween electron and nuclear spins: We have that £1 = 

V7 2 + /? 2 = N = \/E7^F- The o< are 0(1/N), such 
that J2i a f = 0(l/N). When N > 1, then l/N = a < 

= 0(l/y/~N). It is important to notice that in the 
case of an initially fully polarized nuclear state |0) and 
the electron spin in |T) e the system is not affected by the 
exchange interaction terms S-A + and S+A-. 

On the other hand, if the system is initially in the 
state |f) |l)„j the coupling through the Hamiltonian of 
Eq. ([TO]) does not lead to a closed Hilbert space. Then, 
some approximations are necessary in order to obtain so- 
lutions for the evolution of the overall system dynamics. 
Starting from |f 1), we have 



S_i + |Tl) = U)^6«|2y), 



(30) 



where b i3 = a^a^ +a J a i and |2y) = |TT •■■ T< ••■ Tj ••■ ID- 
Defining 



|2>=£ c «l 2 «)' 



(31) 



l<3 



with Cij — bij/N-2, and N2 = we observe that 



i<j 

the state |T 2) is coupled only through 

M- |T 2) = cyoy) |T 1.) = |T 0i) • (32) 

Here, \4>±) can be seen as having one component along 
the state |1) and other component along orthogonal state 
|l p ) , such that 



\<j>i) = (l |<M|l) + vV I0i)-(1|0i) 2 |V, (33) 
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100 200 300 

FIG. 5: Total population Pr of the states 
{|T 1} , |1 2) , |T lp) , \l_2p}}, with N = 10 3 and defect 
distribution a,j — l/VN. 

w here |l p ) = (l/N l p ) (\<j>i) - (1 \fa) |1)), with N lp = 
J (fa \fa) — (1 \4>i) 2 ■ Following this procedure, we have 

S_A+ |T 1 P ) -||2)+ N 2p || 2 p ) ee |T fa) , (34) 

where |2 P ) = jl/N 2p ) (\fa) - (2 \fa) |2)), with \fa) = 

A+ \l p ) and N 2p = \J {fa \fa) ~ (2 \faf- In order to 
implement a semi-analytical description for the overall 
system dynamics we have to truncate the Hilbert space. 
This can be accomplished by approximating the action 
of S+A- over |2 p ) such that 

|| 2 p ) = N 2p |T l p ) + IT Ipp) 

* N 2p \n p ), (35) 

where the state |T lpp) is similarly defined in terms of 
T lp). This approximation can be justified from Fig. [5] 
where we have plotted the total occupation probability 
of states within the subspace of first orthogonal states 
|T 1) , || 2) , |T lp) , || 2 P ). As all of them are eigenstates 
of the J z term, its effect is again a small detuning. These 
results have been obtained by numerical calculation con- 
sidering a 12-dimensional truncated Hilbert space up to 
II 2ppppp)- However, if we use a Hilbert space truncated 
up to || 2pp) , the total occupation probability of the men- 
tioned subspace is practically the same as shown in Fig. [5] 
We have seen before that the typical timescale for the 
exchange of an excitation between the electron and nuclei 
is given by the inverse of the generalized Rabi frequency 

t := n/n = tt/Nq. (36) 

From Fig. [5] we deduce that the dynamics of the com- 
posite system is well described for times beyond r, which 
is important for our subsequent analysis of the quantum 
memory. In particular, for the case N = 10 3 , this time 
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FIG. 6: Evolution of the tangle r e - n (solid line) and popula- 
tion Pj of ||) (dashed line) for a uniformly distributed defect. 
The initial state is |tf) = 1/a/2(|T) + ll» e ® N = 10 3 , 

and r ~ 68t~o. 

is approximately r rj 68to, with To = h/A, and we are 
allowed to restrict the dynamics to the subspace spanned 
by {|T 1) , || 2) , |T lp) • , || 2p)}. This is an important fea- 
ture of our method, since numerical calculations are sig- 
nificantly simplified. 

B. Quantum Information Storage 

In Refs. [13| a long-lived quantum memory based on 
the capability that an electronic spin state can be re- 
versibly mapped into a fully-polarized nuclei state was 
proposed. This dynamics can be expressed as 

(tt|T> + «|i»®|0>-Mi>(u|l)+iu|0>) (37) 

where |1) = (1/Nq) J2i a i This coherent mapping is 
effected by pulsing the applied field from g*fiBB c g 3> Vh 
to B e g ~ for a time r = tt/Nq. Here, we study this 
dynamics in the context of an initial distributed single 
defect in the nuclear spin state. For the transfer process 
to be successful it is necessary that the electronic state 
factorizes from the nuclear spin state, thus requiring that 
at some point of the evolution 

(»m+«\m®m^{}ll-^ti". <38) 

In order to verify if this is possible we analyze the pop- 
ulations of the electron spin and study the entanglement 
between the electron and the nuclear spins. Since the 
overall system is in a pure state we can use as an entangle- 
ment measure the tangle [28j [. tab—^vavb [l — tr (pa)] , 
where pa is the reduced density operator after tracing 
over subsystem B, and z/4 and vb are arbitrary scale 
factors set to va = Vb = 1- 

We first consider the case of a uniformly distributed 
initial nuclear spin excitation. In Fig. [S] we show the 
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FIG. 7: Evolution of the tangle r e _ n for a defect distribu- 
tion peaked at the edge of the quantum dot. Shown are the 
electron spin population P| (dot-dashed line) and the tangle 
(dashed line) for a defect distribution width F — N/50, and 
the tangle for the width F = N (solid line) for N = 10 s and 
the initial state |f) = l/v^dt) + ||)) e ® |l) n . 



electron spin population pi and the tangle between the 
electron and N = 10 3 nuclear spins. We observe that, 
for the case of a uniformly distributed defect, there is no 
time for which the composite state evolves to a factor- 
ized state as in Eq. (|38[) . The electron spin population 
is never maximal or minimal (in both cases the tangle 
would vanish), meaning that there is no high-fidelity in- 
formation transfer to the nuclei. 

We study now the transfer process for the case of a 
localized defect as shown in Eq. . When the defect is 
peaked in the center of the dot, information storage is se- 
riously hindered. However, if the defect is more probably 
located at an extreme of the quantum dot, the coupling 
of the defect-spin to the electron is weak, and separability 
is reached even for rather broad defect distributions, as 
seen in Fig. [71 The necessary conditions for a successful 
transfer are fulfilled after a time t ~ 68to, which is very 
close to the transfer time for an initial fully polarized 
nuclear spin state given by Eq. (136)) . 



The performance of the quantum information storage 
protocol is analyzed now in a more quantitative manner 
by simulating a complete cycle of write-in, storage and 
retrieval. In the first step the electron spin in state 
interacts with the nuclei and after a time r the evolution 
is interrupted and the electron traced out. Then a fresh 
electron in state | J,) interacts with the nuclei and after 
another time r the nuclei are traced out. The fidelity 
F = (ipi\pf\ipi) between this final electron spin state pf 
and the initial state is plotted in Fig. [5] for various defect 
distributions, where F has been calculated by averaging 
on the electron spin Bloch sphere. 



FIG. 8: Fidelity of the quantum memory for N = I0 3 nuclear 
spins and different jo, f < jo < N. 




FIG. 9: Evolution of the population Pf of the electronic spin 
state |t) e for N — 40. The initial state for the electron spin 

is |*> e = 1/^2 (IT) + ll» e and PN = Ef =1 a 3 af |0> <0| <x« 
for the nuclear spins. Solid line: Oj = 1/N; Dashed line: 

a 3 = [c (r/2) 2 /((j - jo) 2 + (r/2) 2 )] 2 , with r = n j = i, 

and C a normalization constant. Dot-dashed line: same as 
for dashed line but jo = N. 



In Fig. [S]we consider the nuclear spins in a mixed state 
where electron spin populations are displayed. This fig- 
ure shows that the electron spin performs very regular 
Rabi oscillations that deviate only slightly from the ideal 
situation with no defect present. In particular the per- 
formance is better than in the corresponding cases of co- 
herently distributed defects. This is to be expected as 
the quantum information storage is sensitive to the co- 
herences of the nuclear spins state, and they are absent in 
this mixed state case. Note that in Fig.[5]the initial num- 
ber of nuclear spin excitations is fixed, at variance with 
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FIG. 10: Evolution ol population Pf of electronic spin 
state |T) e for an initial thermal distribution for each nu- 
clear spin N = 10 and different (k). The initial elec- 
tron spin state is >]/) — l/\/2 + ||}) . For this case, 
(ft) = exp (-PHu) I (1 + exp (-/3hu))) . 



a thermal distribution. In this case the visited Hilbert 
space scales as N 2 and not as 2 N+1 , allowing an exact 
calculation. 

Another possible physical scenario is that of an initial 
nuclear system at finite temperature. Our method is not 
directly applicable to that situation because already the 
initial state covers a large part of the Hilbert space. Nev- 
ertheless, for the sake of completeness, we include a brief 
analysis of this case in the context of a quantum memory 
in the presence of defects. Consider a thermal state for 
each nuclear spin 



Pth 



(39) 



with (3 = I/UbT, where T is the associated temperature. 



The average population of each spin is (k) — 
As for the initial mixed state case of Fig. 



. 1/2 - 
it is not 



possible to define collective states as was the case for pure 
initial states. The results of our direct simulations taking 
into account the whole Hilbert space are shown in Fig.fTOl 
where we show the population of the electronic spin state 
For nuclear polarizations above 95% the contrast of 
the Rabi oscillations of the electron spin deviate from the 
ideal situation by only a few percent. The scaling of the 
error with the particle number N is certainly beyond the 
scope of this article. 



IV. CONCLUSIONS 

We have developed a technique, based on suitable in- 
spection and truncation of the associated Hilbert space, 
that allows to follow the quantum evolution of interact- 
ing systems with inhomogeneous coupling in an accurate 
and controlled manner. As a first general case, we il- 
lustrated the proposed method for the inhomogeneous 
Tavis-Cummings model and showed how the statistical 
properties may change appreciably depending on the spa- 
tial coupling distribution. As a second major example, 
we studied the dynamics of an electron spin interacting 
inhomogeneously with a system of polarized nuclear spins 
in a quantum dot, previously considered as a quantum 
memory candidate. We found that the reliability of the 
storage process depends strongly on the presence and po- 
sition of a single distributed defect in the polarized nuclei 
state. The proposed technique may be easily extended 
to other inhomogeneously interacting systems and may 
prove to be accurate and efficient when interaction times 
are known beforehand, as is the case of pulses. 
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